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ABSTRACT 


The  asymptotic  approximation  method  previously  applied  to  the  case 
of  compressions!  wave  propagation  in  an  inhomogeneous  fluid  medium 
is  carried  through  for  the  case  of  a solid  medium.  Although  the  method 
ia  based  on  the  assumption  of  continuous  variation  of  the  elastic  proper- 
ties, a comparison  of  the  dispersion  curves  computed  by  the  approximate 
theory  with  those  computed  by  the  exact  theory  for  a medium  made  up  of 
two  or  three  homogeneous  layers  indicates  that  the  approximate  theory 
is  fairly  accurate  for  the  normal  modes  of  higher  order  than  the  fifth  or 
sixth.  The  approximation  fails  for  the  modes  of  lowest  order,  which  are 
those  of  greatest  seismological  interest,  but  even  in  this  case  the 
asymptotic  theory,  when  used  in  conjunction  with  the  limiting  forms  of 
the  exact  theory,  has  some  value  for  the  purpose  of  rough  estimation. 
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ASYMPTOTIC  APPROXIMATION  FOR  THE 
ELASTIC  NORMAL  MODES  IN  A STRATIFIED  SOLID  MEDIUM 


1.  INTRODUCTION 

In  a previous  paper,  ^ expressions  were  derived  for  the  dispersion  of  the  normal  modes  of  elastic  waves 
in  a semi-infinite  solid  medium  made  up  of  a finite  number  of  homogeneous  parallel  layers  of  different 
densities  and  elastic  constants.  The  numerical  computation  of  the  phase  velocity  vs  frequency  curves 
from  these  expressions  is  extremely  laborious  if  there  are  more  than  two  layers  and  almost  prohibitively 
so  if  there  are  more  than  three.  In  discussing  the  propagation  of  seismic  surface  waves  in  the  earth’s 
crustal  layers,  particularly  in  the  continents,  it  may  become  necessary  to  employ  more  complex  models  than 
those  that  have  been  treated  previously.  It  is  therefore  of  interest  to  consider  the  afq>licability  to  seismic 
wave  propagation  problems  of  the  asymptotic  approximation  methods  that  have  been  applied  to  fluid  media 
with  continuously  varying  elastic  paramaters.^’^ 

There  is  one  serious  limitation  that  should  be  pointed  out  at  the  outset.  The  approximations  that  we 
shall  use  are  based  on  the  assumption  that  the  properties  of  the  medium  are  slowly  varying  functions  of 
one  coordinate  in  the  sense  that  the  relative  variation  within  a wavelength  is  small.  For  the  frequencies 
of  greatest  seismological  interest  (periods  from  about  0.5  second  to  something  over  a minute),  the  assump- 
tion is  certainly  not  correct  for  the  earth’s  crustal  layers.  Nevertheless  it  may  be  of  some  interest  to  carry 
the  analysis  through  for  the  case  of  a medium  of  slowly  varying  properties  and  then  compare  the  resulting 
normal-mode  dispersion  curves  with  those  obtained  frtm  the  exact  theory  for  discrete  homogeneous  layers. 


2.  SEPARATION  OF  THE  EQUATIONS  OF  MOTION  IN  CYLINDRICAL  COORDINATES 


In  cylindrical  coordinates  the  equations  of  motion  are 


/o3*  u^/ 
pB*  /Bt* 
pB*u^/Bt* 


-dTJ-dr  + (T^~  T^^/r  + dT^g/r-dO  + BT^/B* 
BJ^^/Br  + r„/r  + -dT^g/f-de  + 

+ '^Tgg/'r'de  + Bf^^/Bx, 


(1) 

(2) 

(3) 


Manuscript  received  (or  publication  13  Feb.  1953. 
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where  the  components  o{  the  stress  tensor  ore  given  by 


= kS  + 2/i  3 a,  /Br 

(4) 

^*1  ~ ^ 2/i.3«, /3* 

(5) 

= K6  + 2pCdug/r'bd-  + u,/r) 

(6) 

= /i(3u^/3r  + 3u^3z) 

(7) 

^re  = ~ “5/''  ^ “,A  ^ 

(8) 

+ 3u^/r  3 i9) 

(9) 

and 


S = du^/dr  + u^/r  + ’du^/r'd&  + 3u^/Br  . (jO) 

'fhe  Lamd  elastic  moduli  K and  yu  and  the  density  p are  assumed  to  be  function  of  the  axial  (.aoran  ate 

Z|  only.  Following  Sezawa's^  solution  of  these  equations  for  the  case  of  a homogeneous  medium,  we  start 
with  particular  periodic  solutions  in  the  form 

= Cl)~^  exp  (»Cci<)  cos  n 6 {z)  d J ^ (kr)/dr  ~ {z)  n J ^ (Ar)/r  ] (11) 

= a)"*  exp  (icat)  cos  n 6F^  (z)  (kr)  (12) 

Ug  ~ a;"*  exp  ((OJt)  sin  nd  iF ^{z)  d 1 ^ ih-)/dr  ~ F ^ {z)  n (kr)  /r  ] (13) 

where  n may  be  any  positive  integer  (including  zero),  co  is  the  angular  frequency,  k is  an  arbitrary  num- 
ber which  we  shall  eventually  use  as  a parameter  of  integration,  and  F^,  F^,  and  F^  are  functions  of  z 
to  be  determined.  Since  we  shall  follow  Lamb’s^  method  of  developing  a point  source  solution  by  integra- 
tion of  particular  solutions  of  the  above  form,  the  Bessel  funetioas  of  the  first  kind  (kr),  rather  than  a 
linear  combination  of  functions  of  the  first  and  second  kinds,  are  the  appropriate  radial  functions  in  this 
case.  An  independent  set  of  particular  solutions  may  be  obtained  by  substituting  sin  n 6 for  cos  ra  ^ in 

Eqs.  (11)  and  (12)  and  - cos  nd  foe  sin  nd  in  Eq.  (13).  However,  since  we  shall  not  be  concerned  with 

the  problem  of  representing  particular  types  of  point  sources  by  summation  over  the  azimuthal  characteristic 
number,  n,  the  set  given  hy  Eqs.  (11),  (12),  and  (13)  is  sufficiently  general  for  our  purpose. 
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By  suiMtituting  tii«ae  expreaaions  for  the  dteplacements  in  the  eqaationa  of  motion,  we  obtain  the  follow> 
ing  eqnationa  for  the  functions  F^,  F^,  and  F^ 

F^  + iiF^  + [po;»  - (\  + 2 + {K+  F^  + A =0  d*) 

(\+  2fi)F\  + (\+  2/1)  Fj  + - /iA*]  Fj  - (>.  + /i)  Fj  - \A»  F^  = 0 (15) 

mF,  + /If,  + loco*  -/iA*]F.  =0  (1^^ 

S3  3 

where  the  dot  denotes  differentiation  with  respect  to  s.  If  X,  /x,  and  p were  constants,  the  snbatitntiona 

(17) 

(18) 


^ - Z, 


F,  = A*Z,  - Z, 
2 2 1 


satisfy  Eqs.  (14)  and  (IS)  if  the  functions  Z^,  and  Z^  aatisfy  the  equations 

z\  + [/JcoV(X+  2/x)  - A*]Z^  = 0 (19) 

Zj  + - A*]  Zj  = 0 . (20) 

In  this  case  the  terms  in  Z^,  represent  waves  traveling  with  the  compressional  wave  velocity,  3 - 

[(X  + 2 p)/ p]  and  the  terms  in  Z^  and  F^  (which  satiafies  the  same  differential  equation  as  Z^) 

represent  waves  traveling  with  the  rotational  wave  velocity,  /S  — \.p /p] 

Eqnationa  (17)  through  (20)  are  still  approximately  valid  when  X,  p , and  p are  functions  of  r,  pro* 

• • • 

vided  p/p  A,  X/ XA,  and  p/pk  are  sufficiently  small  to  be  ignored.  This  means  essentially  that  the 
relative  variations  of  the  properties  of  the  medium  within  a wavelength  must  be  very  small.  When  this  is  not 
the  case,  it  will  not,  in  general,  be  possible  to  make  a clear  distinction  between  the  two  types  of  waves. 

In  other  words,  in  an  inhomogeneous  medium  there  will  be  a coupling  between  compressional  and  rotational 
waves  at  every  point  of  the  medium  for  which  p/pk,  etc.,  are  not  negligible. 

Making  the  snbatitntion  of  Eqs.  (17)  and  (18)  and  dropping  the  common  time  factor,  exp  (iox),  the  ele- 
mentary displacements  become 

= o)"*  cos  n 0 [(Zj  - Zj)rf  (kr)/dr  ~ F^n  (Ar)/r  ] (21) 

= w"*  cos  n0  [A*Zj  -Zj]/^(Ar)  (22) 

ug  = o>-»  sin  n 6 [F,  d /,  (Ar)/dr  - (Z,  - Zj)»/,  (Ar)/r]  , (23) 
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and  naing  Elqa.  (19)  and  (20)i  the  atreaa  components  across  a horisontal  plane  become 

cos  - 2/xA*)Zj  + 2 m**  2^]  (24) 

= o)"*  cos  n 6 [{(2  M **  ~ ~ 2 fxZ^  d (!cr)/dr  - pP^n  (*r)/r]  (25) 

= o)"*  sin  b0[m^,  dj^(kr)/dr  - {(2pk^  ~ pco^)  Z^~  2 pZ^nJ  ^{kr)/r]  . (26) 

At  large  values  of  kr,  J ^ (kr)/r  becomes  small  compared  to  ] ^ (kr)  and  d J ^ (kr)/dr,  so  that  the  terms 
in  contribute  appreciablx  only  to  and  the  terms  in  and  Z^  only  to  and  u^.  Thus  at 
large  distances  the  F terms  represent  horisontally  polarized  tranaverse  {SH)  waves,  the  Z,  terms 
represent  transverse  waves  polarized  in  a vertical  plane  (SV),  and  the  Z^  terms  represent  longitudinal 
(P)  waves. 


3.  BOUNDARY  CONDITIONS 

Solutions  representing  a point  source  at  r - 0,  z ~h,  may  be  obtained  by  requiring  the  stress  com- 
ponents to  be  continuous  across  the  plane  z = A,  but  allowing  u^,  and  to  have  discontinuities 
such  that,  when  the  elementary  solutions  are  integrated  with  respect  to  k over  the  interval  zero  to  infinity, 
the  discontinuities  disappear  except  for  a singularity  at  the  source.  The  other  boundary  conditions  are  that 
the  stress  components  T shall  vanish  at  the  free  surface  (z  = 0)  and  at  large  values  of  z 

the  functions  Z^,  Z^  and  must  be  of  ouch  a form  as  to  represent,  with  the  time  factor  exp  (to>t),  only 
waves  propagated  or  exponentially  damped  in  the  z direction. 

Introducing  the  abbreviation  y (z)  = 2 pk^/po>^  - 2{fik/ca)'^,  the  stress  boundary  conditions  at 


the  free  surface  and  at  the  source  are 

[7(6)  - 1]  Zj  (0)  - 7 (0)  Zj  (0)  = 0 (27) 

l7(0)  - 1]  Z,(0)  - 7(0)  Z\  (0)/A»  = 0 (28) 

P,(0)  = 0 (29) 

[y(h)  - 1]  A Z,  - 7(A)  A Z,  = 0 (30) 

(7(A)  - llAZj  - 7(A)  A i, /A*  = 0 (31) 

AP,  = 0 (32) 

where  A signifies  the  discontinuity  in  any  quantity  at  s = A. 
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Let  M^(k,z)  and  (A,*)  be  two  linearly  independent  solutio-  of  E}q.  (19)  which,  with  exp 
represent  waves  propagating  in  the  negative  and  positive  x directions  respectively  for  large  values  of  x. 
Let  U^ikfX)  and  N ^(k,x)  be  the  two  corresponding  solutions  of  Eq.  (20).  The  boundary  conditions  at  the 
free  surface,  the  source,  and  at  infinity  may  then  be  satisfied  by  solutions  of  the  form 


0 < * < A Zj 

= AM^  + 

(33) 

= CM^  + DN^ 

(34) 

= EM^  + FIV^ 

(35) 

h < X Zj 

= 

(36) 

= HN^ 

(37) 

^3 

= //V^  . 

(38) 

By  using  Eqs.  (30)  and  (31)  and  integrating  with  respect  to  k from  zero  to  infinity,  the  discontinuities 


in  the  displacements  at  x - h become 

CD 

A u,  = - o)-*  cos  n ^ / [{A  Z^/y(k)  } d /,  {kr)/dr  + A F,  n /,  (kr)/r]dk  (39) 

0 

00 

A u,  = cos  « 0 / {A  Zj  /r  (A))  /,  (*r)  A»  dk  (40) 

0 

OD 

A sin  « 5 / [ A F,  d /,  (kr)/dr  + { A ZJy  (A)  } n /,  (Ar)/r  ] dk  . (41) 

0 

If  we  set 

AZj  =-o,»y(A)A-iF,  (42) 

AZj  = o.»7(A)A-->S^  (43) 


AF  = 


(44) 
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where  P^,  are  coastants,  the  diacontinuitiea  become 

OD 

A Oj,  = coa  n 6 [P^  (d/dr)  - n Q^/r]  f J ^ (kr)  * dk 

0 

<s 

A = coa  n 6 J I ^ (kr)  A"  * ^ dk 
0 

m 

A Ug  = Bin  a 5 [Q^(d/dr)  - a PJr]  f (kr)  A"*»  dk  . 

0 


The  integral  in  these  expreasions,  conaidered  aa  the  limit  of  / c (Ar)  dk  aa  a approaches 

zero,  vanishes  everywhere  except  at  r - 0,  where  it  haa  an  eaaential  singularity,  thus  satisfying  the 
point  source  requirements. 

From  Eqa.  (33)  to  (38)  we  have 

= (G- B)N^(h)  - AM^(k)  (45) 

AZj  = (H-  D)N^(h)  - CM^(h)  (46) 

APj  = (/  -F)N^(h)  - (A)  . (47) 

By  using  Eqa.  (27)  through  (32)  with  (42)  through  (47),  the  coefficients  in  Z^,  Zj  and  may  be  deter- 
mined as  follows 


A = o^aA-^i  [S,  (7(A)  -l)N^(h)  + P^  y(h)  N^(h)] 

C = -o)2A"*U^-i  [P,  {y(A)-l};V^(A)  + A-»S,  7(A)/V^(A)] 

B = Cj2  p [P^  m (*)  + ;f  (A)] 

D--oo'^  A"^  ‘ p (A)  [P^  y (A)  + Z (A)  ] 

£ = - CO*  A-*i  A‘»  IV^(A) 

P = A-*»  A-i  Af^(O)  Np(k)/Np(Q) 

G = B + [P,r(A)¥„(A)  + {r(A)-l}M^(A)] 

//  = 0 - CO*  A"**  Ay  [P,  (r  (A)  - 1}  (A)  + A-*  S,  r (A)  (A)] 


(48) 

(49) 

(50) 

(51) 

(52) 

(53) 

(54) 

(55) 
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/ = f . (56) 

The  new  eymbole  introduced  in  theae  expressions  are 

~ (57) 

~ ~ (58) 

/ (A)  = [y  (0)  - I]  * (0)  (0)  - [ r (0)/*] » /v„(o)  (o)  (59) 

»'(*)  = 5;*  [y(0)  (r(o)  - iHr(A)  -i}/v^(A)6^  + {r(o)/A}*r(A)  w^(o)/v^(A)/v^(o) 

- iy  (0)  - 1 }2  y (A)  M^(0)  (A)  (0)]  (60) 

^(A)  - A„*  [y(0)  {y(0)  - l}y(A)A"*JV^(A)6^  + {y(0)/AP{y(A)  - D U^iO)  N ^{H)  N^{0) 

-{y(0)-lP{y(A)  - 1}M^(0) /V„  (A)^^(0)]  (61) 

y (k)  = A-i(y(0){y(0) -l)y  (A)  A** /V^(A)  A^  + {y(0)/A}2{y(A)  - l}«^  (0)^^(A)^^(0) 

- {y  (0)  - 1 }» {y  (A)  - 1 } Up  (0)  (A)  (0)]  (62) 

Z(A)  = A-1  A-«[y(0){y(0)  - iKy  (A)  - l) /V^(A)  6^  + {y  (0)/A}*  y (A)  AT^  (0)  (A)  (0) 

- {y  (0)  - 1 }»  y (A)  (0)  Np  (A)  ^^(0)]  . (63) 

The  quantities  and  will,  in  general,  be  functions  of  A,  but  may  be  shown  from  the  differential  Eqs. 

(19)  and  (20)  to  be  independent  of  z.  Since  b^  and  b^  are  the  Wronskians  of  pairs  of  linearly  independent 

functions,  they  cannot  vanish.  Thus  the  coefficients  A,  C and  E ore  finite  for  all  finite  values  of  A, 

while  B,  D,  G and  H have  poles  at  the  zeros  of  f (A)  and  F and  I have  poles  at  the  zeros  of  /V  - (0). 

P 


4.  POINT  SOURCE  SOLUTION  AS  A SUM  OF  NORMAL  MODE  SOLUTIONS 

We  now  intepate  E^s.  (21)*  (22)  and  (23)  with  respect  to  A from  zero  to  infinity  and  evaluate  the  inte- 
pals  in  the  form  of  sums  of  residues  at  the  poles  of  the  integrands  by  applying  Lamb’s*  transfonnation  of 
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the  path  of  integration  in  the  complex  A*plane.  To  simplify  the  resulting  expressions  the  terms  in 
nj^(kr)/r,  which  become  small  at  large  values  of  r,  will  be  dropped.  The  displacements  then  become 

= TTi  cos  n 8 1 [dHl^Uk„r)/dr][df/dk];^  iN (P , Y (kj  + 2 (A„)} 

+ (^J}]  (64) 

u,  = ni  cosnei  r)  [df/dk]  ^k^l^  (z){P^Y{kJ  + S„Z(kJ} 

*m  m '' 

+ /v„  (z)  {P„  r ikj  + s„  ,ir  (AJ}]  (65) 

Ug  = 7T  i sin  « ^ I,  [d  //«)  (r  r)/dr]  {d  N ^Wdk]  (z)  (?„  S^{k)/N^(0)  , (66) 

m 

where  the  Ac^  ’s  are  the  values  of  k at  which  f(k)  - 0 and  the  k^  *s  are  the  values  of  k at  which 
N a (0)  ~ 0.  At  k - k the  functions  W,  X,  T and  Z reduce  to 

P 

lf'(AJ  = [7(0)-  1]  [7(0)  {y(M  - 1}/V^(A)  - 7(M{7(0)  - l}iV^(A)iV^(0)/A^{0)]  (67) 

A(A:„)  = [7(0)  - 1]  [7(0)  7(A)  A (A)  - {7(0)  - 1H7(A)  - 1}A„  (A)/V^(0)/A'„(0)]  (eg) 

T(A„)  = [7(0)  - 1]  [7(0)  7(A)  A-2/V^(A)  - {7(0)  - iH7(A)  - 1> /V^  (A) /V„(0)/A^  (0)]  (^9) 

Z {kj  = [7(0)  - 1]  A-2  [7(0) {7  (A)  - 1}  ^„(A)  - 7(A)  {7(0)  - 1}  (A)  /V^  (0)/A^  (0)]  . (70) 

An  expression  for  the  velocity  potential  in  a fluid  medium  in  a form  analogous  to  Eqs.  (64)*  (65)  and  (66) 
for  the  case  of  a radially  symmetrical  point  source  (n  = 0),  has  been  derived  in  a previous  paper.®  It  was 
shown  that  at  sufficiently  high  frequencies  the  interference  pattern  produced  by  the  superposition  of  a large 
number  of  normal  modes  leads  to  a distribution  of  amplitudes  which  is  the  same  as  that  which  would  be 
computed  on  the  basis  of  ray  geometry.  In  particular,  if  the  variation  of  velocity  with  z is  such  as  to 
produce  a geometrical  shadow  zone,  the  normal  modes  interfere  destructively  in  the  shadow  zone.  At  finite 
frequencies,  the  destructive  interference  is  not  quite  complete  and  there  is  a residual  disturbance  that  is 
given  by  expressions  formally  similar  to  those  occurring  in  the  optical  theories  of  diffraction  at  barriers  and 
caustics.  It  does  not  appear  that  any  new  conclusions  of  much  significance  will  emerge  from  carrying 
through  the  same  analysis  in  the  case  under  consideration.  What  we  do  wish  to  examine  further  is  whether 
or  not  the  asymptotic  solutions  of  Eqs.  (19)  and  (20)  can  be  used  to  draw  any  conclusions  concerning  the 
dispersion  characteristics  of  the  normal  modes  of  low  order,  that  make  the  principal  contribution  to  the  low 
frequency  surface  wave  phases  on  seismograms. 
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Fig-  2-  Velocity  feactioo  for  Caee  IL 
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t 


c < 

A/jj  - {2/nkaJ^  exp  (■*yc,2(*)  + Srri/u] 

(83) 

“i  < c < ««  1 

J. 

s \ 

A/j,  - {2/nka^^  exp  + 5wi/12] 

(84) 

* >«,  J 

a < c < O-l 

1 

(85) 

1 0 1 

/Vj,  - 2{2/vka^  exp  (2  7r»/3)  cos  [k*„  j (*)  " w/4] 

Oj  <r  < aj 

04  < c < aoj 

1 

- 2 (2/wAaj,)'^  exp  Uy„  , (*)  + 5wi/12]  coo  f**a^j(Ol)] 

(86) 

1 < flj  J 

t 

- (2/rrkaJ^  exp  [ -*yo.j  (*)  + 5 7^*712] 

(87) 

* >«J 

c > 

1 

(88) 

, A 

N^  - 2(2/'Tfka^  3xp  (2  7ri/3)  coo  j (x)  “ -n/i) 

* ^ “a  J 

where 

= //l  - Tc/a(«)]^  dx  (89) 

•s 

^<x,  1 (*^  * A ~ [c/a  U)1  * dt  (^^ 

Z 

*2  

*0,2  (*^  ■ ^ /Tc7oT(*)P^  rfr  (91) 

t 

and  for  c < <Xj  we  take  o^  " ®j  ~ 

To  the  same  order  of  approximation  the  derivative  is  given  by 
Case  I 


c < a 

0 

/Vj,  -*  (2  ka^/v^  exp  t -ky^(t)  ~ 7ri/12] 

(92) 

“’“.l 

. 1 

- (2  ka^TT^  exp  [ -ky^{z)  - 77iA2] 

(93) 

X > 0 J 

/Vjj  -•  2 (2  ka^Tt^  exp  (2  7ri/3)  sin  (x)  - 77/4]  . 

(94) 

X < 0 J 
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Case  n 


1 


c < 

/V^  - (2  exp  [ “ ky^  ,,  (*)  “tti  /12] 

(95) 

a,  < c < (to"! 

4_ 

^0  - *■«  “P  ^ * ^o.a  ■ ^*’/12i 

(96) 

* >a,  J 

a,  < c < 

1 0 

. 1 

/Vj,  - exp  (2  7Tt/3)sin  [kx^^^ o (r)  - tt/4] 

(97) 

0,  < X < Oj 

ttj  < C < dA 

*<«l  J 

J 

yv^  - 2 (2  kn^r^*  exp  t k y^  j (*)  - tti  /12]  coe  ikx^^  , (Oj)J 

(98) 

- (2  kaJrtS^  exp  [“kyj,^  o («)  - rr»  /12  ] 

(99) 

* >«J 

/Vj,  - 2 (2  kr^rr)*  exp  (2Tr«/3)sin  [kx^^  ^(z) -7r/4]  . 

(100) 

* <a,J 

He  foBctions  N ^ sad  are  given  by  the  aame  expceasiona  witb  fi  subatitnted  for  a.  Correspond* 
ing  values  of  r at  which  = e are  denoted  by  b.  The  characteristic  Eq.  (72)  for  the  normsl  modes  of 
Rayleigh  type  then  takes  the  following  forms. 

Case  I 

[y(0)-l]»  - [y(0)]  * /[l  -(c/a^)*)  [i  -(c/y0j)*]  = o UOl) 

/So  < c < Oo  [y  (0)  - 1]  » + [y(0)j  » /[l  - (c/tto)*]  [(c//So)»  - 1]  tan  tk,  (0)  " tt/4]  = 0 (102) 

«0  ^ « (y  (0)  - 1]  * - [y  (0)]  * /[(c/Oo)*-i][(c//0o)*-i] 

tan  *„(0)  " w/4]  tan  [k^  »^(0)  ”7r/4]  =0  . (103) 

Case  n 

c < [y  (0)  -1]  * - [y(0)3  * /[l  -(c/Oo)*]  tl  -(c//0o)*]  = 0 (104) 

y8l<e</8o<«i  |[y(0)-j]  * - [y(0)]  V[l-(c/ao)»]  [1  -(c//5o)>]|cos  [k^  s^  ,(kj)J  =0  (105) 

/So  <C  < «1  Cr(0)  -1]  * + [y(0)l » /[l-(c/ao)»]  [(c//3o)*  - l] 

ton  f*»  */j.,(0) -w/4]  =0 


(106) 
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Oj  < c < 


®o  ^ [y(0)]  * / [ 1 [(c//3,j)*  - l] 


tan 


(«i)]  =0 


(107) 


c > a. 


[y(0)-l]»-  ty(0)]  V [(c/a/  - l]  [(c/;3/  - l]  tan  ,(0)  -w/4] 

tan  [*^*^  j (0)  - TT/4]  = 0 . (108) 


Equations  (101)>  (104)t  and  the  first  factor  of  Eq.(lOS)  are  the  same  as  the  equation  for  the  velocity  of 
Rayleigh  waves  on  a homogeneous  medium  having  velocities  a = a^  and  /i  ~ This  equation  deter- 
mines a single  real  value  of  the  phase  velocity,  c,  for  all  values  of  k.  Actually,  of  course,  this  value  of 
c is  merely  the  asymptote  approached  by  the  Rayleigh  wave  phase  velocity  at  high  frequencies,  so  that, 
the  present  approximation  is  inadequate  to  deal  with  the  dispersion  of  the  mode  of  lowest  order* 

The  characteristic  Eq.  (73)  for  the  normal  modes  of  SH  type  has  the  following  asymptotic  fonns. 


Case  I 

c < /3.  No  roots  exist 

0 

c > sin  x^(C)  “ TT/4J  = 0 

or  ~ 1/4)  « = 0,1, 2, 3 (109) 

Casen 

c < /3^  No  roots  exist 

< c </3^  cos  [*;*^,(frj)]  =0 

or  ^ 1/2)  m = 0, 1,2,3 (110) 

< c sin  [*„|*y},j(0)  ” t'/A]  = 0 

* *«*^,(0)  = w(m  + 1/4)  m = (Ill) 


6.  APPUCATION  OF  APPROXIMATC  THEORY  TO  A LAYERED  MEDIUM 

Hie  development  of  the  approximate  ezpraasions  for  the  fnnctiona  and  N ^ bvolvea  the  aasnmption 
that  the  velocities  and  their  first  derivatives  are  eostinnoss  fnnctieos  of  s so  that,  in  general,  the  present 
theory  is  not  qiplicable  to  a medinm  made  np  of  discrete  layers.  If  the  diseostinnities  in  the  vsloeitios  and 
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Fig.  3.  Comparfaiaa  of  diapersion  corvea  for  Rayleigii  Modea  in  a two  layered  nedinin. 


the  density  are  snail,  however,  the  anplitndea  of  the  waves  reflected  from  the  boundaries  will  be  small  ex- 
cept for  angles  of  incidence  greater  than  the  critical  angle.  If  we  can  ignore  the  boundary  reflections  at 
angles  less  than  the  critical  angle,  there  is  no  fundamental  distinction  between  the  cases  of  continuous  and 
discontinuous  velocity  variation. 

In  the  case  of  monotoaically  increasing  velocity,  the  quantities  (0)  and  x^(0)  become  for  discrete 
layers,  each  of  constant  velocity. 


*a(0) 


<■1 


*^(o) 


fi.<c 


<-l 


(112) 


- 1 


(113) 


a 
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Fig.  4-  Comparison  of  dispersion  curves  for  SH  modes  in  a two  layered  medium. 


where  d.  is  the  thickness  of  the  tth  layer  and  the  summations  are  extended  over  all  layers  for  which 
or  respectively  are  less  than  c.  If  the  velocity  of  successive  layers  decreases  and  then  increases, 

we  find  the  same  Eqs.  (112)  and  (113)  for  ^(0)  and  x^.  ^(0)  when  c > and  c > respectively. 
When  c OL^  or  c < the  summations  extend  only  over  the  layers  for  which  < c or  < c. 

Phase  velocity  dispersion  curves  have  been  computed  according  to  the  present  approximate  theory  and 
according  to  the  exact  theory  for  a two*layered  medium  (the  second  being  senMnfinite)  having  the  following 
properties: 


= 1.81 

= 1.37 

aj//3,  = 2.44 

= 1.11  ■ 

The  results  are  plotted  in  dimensionless  form,  c/yS^  versos  kd,  for  the  first  seven  normal  modes  of 
Rayleif^  type  in  Fig.  3,  and  for  the  normal  modea  of  SH  type  in  Fig.  4.  The  exact  and  approximate  curves 


I 


Fif.  S*  CoMparisoa  of  dioporaiaa  corveo  for  SH  medoo  to  • ttooe  lajrored  modioa,  Case  I. 

•n  coaqwrad  for  tke  Sff  aodM  !■  two  S*UifM«<l  aedio  fai  Figs.  5 ud  6.  Tliooe  shown  in  Fig.  5 are  for 

ifcoeoM 


= 1.19 

- 1.00  = 1.56  , 

/8,//8j  = 1.37 

= 1.11 

aad  thooa  shown  in  Fig.  6 illnatrsts  dio 

offset  of 

a low  eslocitjr  socoad  lajrer 

y8,//3,  = .988 

Pjf\ 

- 1.00  d^dj  = .871 

fijfi.  = 1.87  Al/Al 


1.11 


i 
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Fig.  6*  Comparison  of  dispersion  curves  for  SH  modes  in  a three  layered  medium,  Case  II. 


From  these  comparisons,  it  appears  that  the  approximate  theory  gives  a fairly  good  representation  of  the 
phase  velocity  dispersion  curves  for  the  normal  modes  of  order  higher  than  the  5th  or  6th,  except  in  the 
immediate  neighborhood  of  the  low  frequency  cut  offs.  For  the  lower  modes  the  approximation  is  not  good 
enough  for  quantitative  purposes.  In  particular,  in  studies  of  seismic  surface  waves  one  will  usually  want 
to  compare  an  observed  variation  of  gyoup  velocity,  U ~ d{ck)/dk,  with  period,  with  the  theoretical  gyoup 
velocity  curves  for  the  modes  of  order  m = 0 or  I computed  for  some  assumed  velocity  stratification.  It 
is  quite  evident  that  the  approximate  curves  for  the  cases  illnstrated  here  would  not  be  adequate  for  this 
purpose. 

It  is  believed,  however,  that  the  approximate  theory  will  be  found  to  have  some  value  for  the  purpose  of 
qualitative  estimation  even  for  the  modes  of  lowest  order.  For  example,  in  the  case  shown  in  Fig.  3,  we 
know  that  the 'correct  curve  for  m ~ 0 should  approach  the  Rayleigh  velocity  of  the  semi-infinite  layer  os 
k approaches  zero,  and  from  the  asymptotic  form  of  the  exact  theory  for  small  valhes  of  k,  we  have  an 
expression  from  which  the  behavior  of  this  carve  lor  kd  <<  1 can  be  computed  without  difficulty.  We  also 

i 


tt 


know  that  for  kd>  10  the  phase  velocity  of  this  mode  will  be  very  close  to  the  Rayleigh  velocity  of  the 
upper  layer.  With  theae  limitiag  valnea  fixed  and  knowing  that  the  approximate  curve  approaches  the 
coirect  curve  for  some  intermediate  value  of  c,  the  general  form  of  the  curve  can  be  estimated  within  a 
moderately  broad  band  of  uncertainty.  In  the  case  under  consideration,  we  could  have  inferred  by  auch  a 
process  that  the  minimum  gyoup  velocity  would  probably  occur  at  some  value  of  kd  between  2 and  4.  The 
actual  minimum,  obtained  by  a gyaphica)  method  from  the  plotted  “exact”  curve,  is  at  kd  = 3.05.  In  more 
complex  cases,  where  the  labor  of  computing  the  dispersion  curves  by  the  exact  theory  may  become  prohibi- 
tive, even  so  loose  an  estimate  as  this  may  be  useful  in  limiting  the  range  of  possible  structural  models  that 
adght  be  used  to  explain  a given  set  of  obaerved  dispersion  data. 
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